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An expanding closed universe filled with radiation can either recollapse or tunnel 
to the regime of unbounded expansion, if the cosmological constant is nonzero. We 
re-examine the question of particle creation during tunneling, with the purpose of 
resolving a long-standing controversy. Using a perturbative superspace model with a 
conformally coupled massless scalar field, which is known to give no particle produc- 
tion, we explicitly show that the breakdown of the semiclassical approximation and 
the "catastrophic particle production" claimed earlier in the literature are due to an 
inappropriate choice of the initial quantum state prior to the tunneling. 

I. INTRODUCTION 

The tunneling approach to quantum cosmology proposes to view the creation of the 
Universe as a quantum tunneling event @, |3|- A semiclassical picture based on the 
Wheeler-DeWitt equation describes tunneling from a state of vanishing size ("tunneling 
from nothing") to a closed inflating universe. Alternative proposals have been put forward 
by Hartle and Hawking |Q and Linde ||. For a recent critical review see, e.g., 0. In the 
present paper, we shall focus on the tunneling picture. 

The process of tunneling from nothing can be thought of as a limit of tunneling from 
a closed recollapsing universe of very small but nonzero size ("tunneling from something"). 
Consider a closed FRW universe filled with a vacuum of constant energy density and some 
matter or radiation. Classically, such a universe has two possible types of evolution. It can 
expand to a maximum radius and recollapse, or it can contract from an infinite size, bounce 
at a minimum radius and then reexpand. But in quantum cosmology there is yet another 
possibility. Instead of recollapsing, the universe can tunnel through a potential barrier to the 
regime of unbounded expansion. The semiclassical tunneling probability does not vanish in 
the limit when the maximum size of the initial universe shrinks to zero. The corresponding 
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wave function describes quantum nucleation of an inflating universe from nothing. 

The question of the particle content in the universe created by quantum tunneling has 
been studied in the literature by including the inhomogeneous modes of quantum fields in 
the wave function. Conflicting claims of excessive particle production during tunneling [Q, [7J, 
on the one hand, and of essentially no particle content in the inflating universe || ||, [TO], 
on the other hand, have been advanced. In this paper we intend to resolve this controversy 
and explain the origin of the conflicting results. 

Rubakov (2| (see also (?]]) considered scalar particle production in the process of "tun- 
neling from something". He found that, for a generic initial state of the universe, the wave 
function is significantly affected by the apparent growth of the excitations of the scalar field 
during tunneling. Rubakov interpreted this result as an excessive particle production during 
tunneling and a breakdown of the semiclassical evolution. However, we shall argue that this 
result is a consequence of an inappropriate choice of the quantum state of the universe and 
that in fact there is only a finite, if any, particle production during tunneling. 

Our main points can be summarized as follows. An arbitrary initial state of the recollaps- 
ing universe can be interpreted as a superposition of semiclassical geometries with certain 
amplitudes, each branch being a classical universe with a certain quantum state of the scalar 
field. In a generic superposition, essentially all excited states of the scalar field will be repre- 
sented, perhaps with small amplitudes. In this case, some of the semiclassical branches will 
have such a large energy due to the excitations of the scalar field that no tunneling will take 
place: in these branches, the initial universe will not recollapse but will continue expanding. 
Non-tunneling branches give a potentially larger contribution to the wave function at large 
scale factors because they are not exponentially suppressed, compared with the tunneling 
branches. If the amplitudes of such high-energy branches in a given superposition are suf- 
ficiently large, they will give a dominant contribution to the wave function of the inflating 
universe, while the contribution of the tunneling branches will be negligible. 

If one considers the tunneling from a recollapsing universe of nonzero size ("tunneling 
from something"), then one is free to choose the initial quantum state of that universe. 
We shall show that a consistent semiclassical picture of the tunneling universe can only be 
obtained with an appropriate choice of the initial quantum state of the recollapsing universe. 
This state should be such that both the recollapsing and the inflating universe belong to the 
same semiclassical branch of the wave function. Rubakov et at. 0, [7]] have chosen a quantum 
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state that does not satisfy this condition. On the other hand, Refs. || § have considered 
"tunneling from nothing" (the limit of zero size of the recollapsing universe) using a quantum 
state that describes a single underlying semiclassical geometry. In the present paper we shall 
explicitly construct such quantum states for the case of "tunneling from something". The 
difference in the choice of the quantum state is the first cause of the discrepancy in the cited 
papers. 

The second problem with the results of Rubakov et al. is their particle interpretation of the 
wave function. A particle interpretation of a quantum field theory requires a fixed, classical 
background metric. With a generic choice of the quantum state of the universe, the Wheeler- 
DeWitt wave function will describe a quantum superposition of different geometries, rather 
than a single classical geometry. In this case, the interpretation of the Wheeler-DeWitt 
wave function that Rubakov et al. used to obtain the particle content (the formalism of the 
instantaneous Hamiltonian diagonalization) is not justified. If the particle content is inferred 
from the wave function as if there exists a unique underlying semiclassical spacetime, then 
one is lead to erroneous conclusions about the breakdown of the tunneling process and about 
the excessive particle production. 

The paper is organized as follows. In Sec. 2, we review the "perturbative superspace" 
approach to the Wheeler-DeWitt equation. We demonstrate that a Gaussian ansatz for 
the WKB wave function employed in Refs. M, |!| corresponds to an "instantaneous squeezed 
state" in the formalism of Rubakov et al. In this sense we find a formal agreement between 
these calculations. In Sec. 3, we consider a massless conformally coupled scalar field, in 
which case the Wheeler-DeWitt equation is separable and it is well known that there is no 
particle production ||1 1|| . Assuming a squeezed initial quantum state for the scalar field, we 
show by an explicit calculation that the wave function becomes dominated by high-energy 
states far enough under the barrier, and that the interpretation of Rubakov et al. would 
indicate a "catastrophic particle production" during tunneling. Details of the calculation 
are given in Appendix A. In Sec. 4, we discuss what we believe to be the correct physical 
interpretation of the results. The more complicated case of a massive scalar field will be 



presented in the companion paper 12 



4 



II. SEMICLASSICAL PERTURB ATIVE SUPERSPACE 

We consider a homogeneous (closed) FRW universe with a conformally coupled scalar 
field <p. The metric is homogeneous on 3-spheres, 

ds 2 = dt 2 - a 2 (t) dQ 2 3 , (1) 

where a (t) is the scale factor. The scalar field is not homogeneous and may be expanded 
in 3-spherical harmonics, 

*fo*) = ^£W*W*)- ( 2 ) 

n,l,p 

Below, only the index n = 1, 2, ... will enter the equations, and we shall suppress the indices 
I, p of the modes Xnipif). The summation over degenerate indices I, p spans I = 0, 
n — 1 and p = —I, / and introduces an extra factor n 2 which we shall insert in explicit 
calculations below. 

The model is described by the classical action 

d\^~g { -Lr + I (d^f -^H 2 -^-±-rA. (3) 



16tt 2 v 8tt 2 12 

Here, the parameter if represents the vacuum energy (the cosmological constant), R is 
the scalar curvature, m is the mass of the scalar field, and we are using Planck units, 
G = h = c = 1. In the Schrodinger picture of the "perturbative superspace" approach 
to quantum gravity 0], the wave function of the universe ^ (a, {%„}) after appropriate 
rescalings of the parameters (see [§] for more details) satisfies the Wheeler-DeWitt equation 

* (a, {Xn}) = 0. (4) 



^-V(a)-][> 



da 2 

Here, the scalar field Hamiltonian for the n-th mode is 

H n = h 2 ^ - (n 2 + m 2 a 2 ) xl (5) 

dxi 

and 

V (a) ee CL 2 — H 2 a\ (6) 

We have written out the Planck constant h in Eq. (|) to make the WKB approximation 
more explicit below. [The factor h is merely a formal bookkeeping parameter since H = 1 
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in Planck units.] The Hamiltonian of Eq. (|5|) describes a harmonic oscillator with an a- 
dependent frequency 

u n = UJ n (a) = Vn 2 + m 2 a 2 . (7) 

As in Ref. 0, we shall include in addition to the scalar field a homogeneous radiation 
component with energy density 

Pr = a~ A e r , (8) 
where e r is a constant parameter. This amounts to replacing 

H 2 - H 2 + % (9) 

Qi 

and therefore Eq. (|J) still holds with 

V(a) = a 2 - H 2 a 4 - e r . (10) 
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Figure 1: "Tunneling from something": creation of the inflating universe by tunneling. 

The physical picture of the universe in this model is illustrated in Fig. |1|. The Wheeler- 
DeWitt equation in the (a, Xn) space is formally similar to a stationary Schrodinger equation 
for a quantum-mechanical particle in a two-dimensional potential. A small closed universe 
filled with radiation of energy density e r a~ 4 is expanding and recollapsing in the neigh- 
borhood of a = 0, and an inflating universe with scale factor a 2 is created by tunneling 
through the potential barrier. The maximum scale factor a\ of the recollapsing universe is 
determined by the energy parameter e r . 

A typical solution of the Wheeler-DeWitt equation that represents the creation of an 
inflating universe by tunneling is qualitatively analogous to the quantum-mechanical tun- 
neling wave function with the above potential. In the regions outside the barrier, the wave 
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function is oscillatory (in the a direction) and generally has the form of a linear combination 
of incoming and outgoing waves, 

* (a, {Xn}) = h ({Xn}) e 45l(a) + h {{Xn}) e-^ a \ (11) 

The tunneling boundary condition specifies that the wave function should contain only the 
outgoing wave in the domain a > a 2 - In this domain, we can interpret the wave function 
as describing a classical spacetime if we introduce a semiclassical time variable which is a 
function of a III, [TH|. A customary ansatz is to introduce the conformal time by 



dr= da . (12) 



In the underbarrier region, a,\ < a < a 2 , the wave function is a linear combination 
of growing and decaying real exponentials; the conformal time becomes imaginary. The 
underbarrier region a% < a < a 2 does not correspond to a classical spacetime. 



A. The Gaussian approximation 



To obtain an approximate solution, one can substitute the following WKB-motivated 
Gaussian ansatz into Eq. (13), 



* (a, {Xn}) = exp 



-h S{a) -2H 



(13) 



and neglect the terms of order 0(xt)- (This ansatz was introduced in Ref. (T^| in the 
context of tunneling in two dimensions.) Equation (|i~3|) approximates the wave function 
only near Xn — 0, where higher powers of Xn are negligible. We shall refer to the wave 
function of Eq. (|i~3"D simply as the Gaussian solution. In the companion paper [|TJ| we shall 
give a detailed analysis of the applicability of this approximation. 

The functions S (a), S n (a) in the ansatz are to be found from the equations 

{S') 2 -V(a) — hS" + h Sn = 0, 



(14) 

S'S' n -Sl + n> + m>a> - ^ = 0. (15) 
The WKB approximation consists of neglecting the terms of order 0(h) in Eqs. (|T4])-(|To|). 



[Note that the last term in Eq. (|i~4D can be interpreted as the backreaction of the scalar field 
excitations on the background metric] 
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In the underbarrier region a\ < a < 02, we introduce the Euclidean conformal time 
variable r by 

, , f a da , , 

Tia)s Lvm- ,16) 

[A straightforward analytic continuation of this definition is to be used in the classically 
allowed regions.] Then we obtain the equations 

S(a) = ± J da^V (a), (17) 
d -^ = Sl-ul (18) 

A perhaps unexpected property of the Gaussian solution is that the real part of the 
function S n (a) obtained from Eq. (|18| ) may become negative at some values of a under the 
barrier, making the Gaussian wave function grow at large Xn- Whether or not this happens 
is determined by the initial value S n (ai) at the first turning point. In previous work Q, the 
condition 

ReS n >0 (19) 

was motivated by the requirement that the wave function be finite at large Xn- This condition 
has been used in Refs. || |10| and yielded physically reasonable results. However, it should 
be noted that, strictly speaking, the occurrence of Re S n < is not necessarily problematic, 
since the Gaussian approximation should only be applicable at small Xn- We shall see 
below that the onset of the "catastrophic particle production" as claimed by Rubakov et 
al. is directly related to the change of sign of ReS n (a). Our considerations will clarify the 
physical interpretation of this phenomenon and justify the condition of Eq. (|TTf). 



B. The instantaneous diagonalization picture 

Rubakov et al. 0, 0] have used the method of instantaneous Hamiltonian diagonalization 
to solve Eq. (||). This is equivalent to expanding the wave function \1/ (a, {Xn}) in the eigen- 
states (a, Xn) of the a-dependent Hamiltonian of Eq. (|[). The normalized instantaneous 
eigenstates can be taken as 

^^-SSWW^M-^)- (2o) 
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where Hk (x) are Hermite polynomials and uj n (a) is given by Eq. ([<]). The Xn-dependent 
part of the wave function is decomposed into a superposition of instantaneous excited states, 

oo 

*n (a, X n) = ^ («) 4 n) (a, Xn) , (21) 

A:=0 

with unknown a-dependent coefficients cj^ (a). The complete wave function is a product 
of these factors, 

^(a,{ X n}) = e- s ^l[^ n (a,Xnj 

n 



An approximate solution obtained in Ref. [§] is of the form 

(a) = [C» (a)]* , (22) 

where ( n (a) is a certain explicitly obtained function, which quickly grows with a when 
starting at ( n (ai) = 0, making the amplitudes (a) large. However, the expansion 
of Eq. ( f2ll ) with the coefficients from Eq. (^) is meaningful only when |£ n (a)| < 1. In 
Ref. [§], the coefficients are interpreted as amplitudes for the mode Xn to be in k-th 
excited state, with the conclusion that such large amplitudes demonstrate a catastrophic 
particle production and a breakdown of the perturbative and/or the WKB approximation. 
Although the coefficients C^f 1 should not be interpreted as real particle numbers in the 
tunneling regime, it is clear that the backreaction cannot be neglected if become large. 

The choice of the initial value ( n (ai) = was motivated in Ref. by the intention 
to investigate the particle production during tunneling and to start with very few or no 
particles in the recollapsing universe. Thus the state with no particles instantaneously at 
a = a±, i.e. C (ai) = 1, (ai) = for k > 1, was chosen. However, the definitions of 
particles and of the vacuum in an expanding universe are notoriously ambiguous (see, e.g., 
|T7|| ), and the method of instantaneous diagonalization is known to give unphysical results 
in some cases (11, An unambiguous definition of particles is possible only in certain very 



special models. One such model is a massless conformally coupled scalar field, for which 



there is strictly no particle production ||1 1|| . We shall show in the companion paper [12[] that 
the state chosen by Rubakov et al. cannot be considered a vacuum state. In the next section 



we shall interpret the states described by Eq. (|22 ) as squeezed states. 
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C. Instantaneously squeezed states 

A squeezed vacuum state |C) of a harmonic oscillator can be defined using the creation 
operator a) and the vacuum state |0) as 



|C) = (l-|C 2 |) 3 exp 



|0> • (23) 



Here C is the "squeezing parameter", a complex number satisfying |£| < 1. The state is 
normalized with the given prefactor, so that (CIC) = 1- Since 

(J) k \0) = Vk\\k), (24) 

we find that the squeezed state is the following superposition of the excited states, 

lO^l-ICI^E ^ftPV ). (25) 

fc=0 

[The normalized wave functions for the excited states are given by Eq. ( p0| ) and the choice 
of phases in that equation is consistent with Eq. (p4|).] 

Suppose that at a particular moment of time the wave function of the harmonic oscillator 
with frequency u n is a Gaussian, 

4>G (Xn) oc exp ( --^xl J , (26) 



2 

and assume that ReS n > 0. We can represent the wave function of Eq. (|26() by a superpo- 
sition of excited states, 

oo 

\^ G ) = Y,C k \k) . (27) 

fc=0 

Only even k = 2p will have nonzero amplitudes C*.. Using the following formula for the 
integral of a Hermite polynomial (transformed from Eq. 22.13.17 of |[T9| ), 

p+l/2 



H 2p (x) exp [~ bj pj dx = V^ (2p - 1)!! (2 - bf 



(28) 



where b > 0, we obtain, up to a p-independent factor, 



^-SA ' (2p-l)M 

C2p K l^T^J VST' 



Comparing this with Eq. (E^) and using the identity 



(^" 1 ) !! = ^ (30) 
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we find that the a-dependent decomposition of the Gaussian wave function of Eq. (|26|) into 
instantaneous excited states is exactly the same as that of a squeezed vacuum state with 
the (a-dependent) squeezing parameter 

S n (a) + u n (a) 

Incidentally, from Eq. ([15]) it follows that the function ( n (a) in the underbarrier region 
satisfies the equation 

f = (32) 

This equation will be useful in the companion paper ||12|| . 

Thus we have identified the function Q n (a) from Eq. fl22|) as the instantaneous squeezing 
parameter. This allows us to relate the change of the signature of the Gaussian solution with 
Rubakov's "catastrophic particle production". The instantaneous mean occupation number 



in a squeezed state of Eq. fl25|) is 

(N) = (1 _ \ C 2\^ IC|2fc (2p)! 2 /: = IC|2 (33) 

The mean occupation number (in one mode) becomes formally infinite when \(\ grows above 
1. However, Eq. (|3l|) can give |Cn(o)| > 1 only when ReS n (a) < 0. Therefore, the mean 
occupation number remains finite, and an apparent catastrophe is avoided, as long as Eq. (|19D 
holds. 

We still have not given a physical motivation for the condition of Eq. ( |19| ) . In the next 
section we shall address this issue by considering an exactly solvable example. 

III. EXPLICIT RESULTS FOR A MASSLESS FIELD 

In this section we consider the model of a tunneling universe with a conformally coupled 
massless scalar field. It is well known that there is no particle production in this model 



Till . However, the arguments of Rubakov et al. would still suggest a "catastrophic particle 
production" for some states of the field. To understand the origin of this discrepancy, it 
is instructive to compare the Gaussian solution with the exact wave function obtained by 
separation of variables in the Wheeler-DeWitt equation. 
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A. Solution by separation of variables 



In the massless case, the Wheeler-DeWitt equation is 

2 



If 



da? 



n 



*(a,{ X n}) =0. 



(34) 



[Here we have inserted the degeneracy factor n 2 .] The variables in Eq. fl3"4p separate. The 
separable solutions are of the form 



(35) 



A general wave function is a linear combination of such solutions. 

A single separable solution of the form ( |35| ) may be parametrized by integer occupation 
numbers p n > 0, indicating the excitation levels p n of the modes Xn- The eigenfunctions in 
Xn are given by Eq. (|20|) with the substitutions uj n = n and k = p n . For a given set of the 
occupation numbers {p n }, the eigenfunction in a satisfies 







( fl ) = 0. 



(36) 



Note that since the number of modes is infinite, the above sum diverges even if all p n = 0; 
this is the divergence of the zero-point vacuum energy. To obtain a meaningful solution, 
we assume that the zero-point energy divergence is absorbed into the radiation density 
parameter e r , and that only finitely many of p n are nonzero. Then the equation for the 
eigenfunction ip (a) becomes 

d 2 



h 2 —-V(a) + hY2n- 
da 1 



Pn 



( a ) 



0. 



(37) 



The 2n 3 p n term represents the backreaction of the scalar field excitations on the background 
geometry. This term becomes significant if we consider excited states with large p n . 
Equation (|3~7D can be solved using the WKB approximation, 

S(a) 



^{ P u} (a) = exp 



h 



(38) 



with S (a) satisfying the equation 

dS 



— = ±JV(a)-hJ2^ 3 Pn- 



(39) 
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We have kept the O (H) term in Eq. (|39|) because we would like to allow arbitrary combina- 
tions of excitation levels, which may lead to large values of ^ n 2n 3 p n . 

Because of the two possible signs at the square root, we obtain two branches corresponding 
to the growing and the decaying solutions under the barrier, 



S 9 



f 1 lv{a)-hJ2^ 3 Pn da, S d = H lv (a) - n^2n 3 p n 



da 



(40) 



Here, the boundaries a± t 2 of the classically forbidden region a± < a < a 2 are the two positive 
roots of the equation 



Via) 



Pr, 



0. 



(41) 



The boundary conditions for the a-dependent part of the wave function are a normaliza- 
tion condition, e.g. ip (a = 0) = 1, and the tunneling boundary condition at large a. 

In the case of zero occupation numbers (the vacuum state), the general (WKB) solution 
of Eq. fl34|) under the barrier is a linear combination of the two branches, 



* (a, Xr. 



C„e « +C d e 



exp 



1 

2h 



n 



(42) 



The subscript in Sq 9 signifies that they are the vacuum solutions given by Eq. (|l]) with all 
p n = 0. The relation between the constants C g and Cd is to be obtained from the boundary 
condition at a = oo and will not be important for what follows. 
An excited state of the mode Xn' is given by the wave function 



* (a, Xn) oc H. 




(P) „ (P) 

ft + C d e ft 



(43) 



where the functions S 9 ^ are given by Eq. ( |40| ) with all the occupation numbers p n equal to 
zero except p n i = p. These functions are different from the vacuum functions because 
of the backreaction of the scalar field excitations on the metric. 

Equations (|i"2"D-(fi"3"l) apply in the under-barrier region, a\ < a < a 2 ; analogous expressions 
can be written for the other regions. 

Because of the separation of variables, the definition of the vacuum and the particle 
interpretation of the wave function are unambiguous. A given wave function \1/ (a, {x n }) 
is decomposed into separable solutions of the form of Eq. fl35|). Each branch corresponds 
to a semiclassical universe with a fixed set of occupation numbers {p n } and a modified 
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background geometry. Since the background geometry is affected by the excitations, a linear 
superposition of several such semiclassical wave functions with sufficiently different sets of 
occupation numbers {p n } may not correspond to a universe with a definite semiclassical 
geometry or particle numbers. However, any classical observers in such a universe will find 
themselves in certain semiclassical branches where the spacetime is fixed and the particle 
numbers remain constant. In this sense, there is no observable particle production in the 
massless model, for any choice of the quantum state of the universe. 



B. Solution in the Gaussian approximation 

If So (a) satisfies Eq. (|39D with either of the signs at the square root, then Eq. flilf ) can 
be solved in quadratures (for any V (a)). We again introduce the conformal time variable r 
by Eq. (^). The general solution of Eq. (0) with the identification u n = n (as appropriate 
for the massless case) is 

e -2nr{a) _|_ _g 

Sn(a) = n ^ 2nT{a) _ B , (44) 

where B is an arbitrary constant of integration. We can rewrite the solutions for the growing 
and the decaying branches as 

Qg( l + B g e 2nT ^ ^ l + B d e- 2nT ^ 

S 9 ri (a)=n — - — 7— r- , S„(a)=n „ , , , (45) 

where we have introduced the constants B g and B d for these two branches. (We suppressed 
the index n at the constants B g , B d for brevity.) 

The choice of the constants B g ^ is determined by the initial condition S n (ax) and by 
matching of the branches at the second turning point 02. It is clear from Eq. fl4"5[) that 
some choices of the constants will lead to negative values of S n (a) under the barrier. A 
special choice is B g = B d = 0; in this case, the function S n (a) is constant, and the vacuum 



solution of Eq. (|4"2 ) is recovered. For other values of B, the Gaussian solution corresponds 
to a squeezed state. The corresponding squeezing parameter at a — a\ is ( n (ai) = B, as 
follows from Eq. fl3"T|) with the identification uj n = n appropriate for the massless case. 
The squeezing parameter ( n (a) satisfies Eq. (|32|) where uj n = const, 

f^2<„. (46, 
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For a given initial squeezing parameter ( n (ai), we find that, at a point a = a* under the 
barrier, |£ n (a*)| = 1 and therefore ReS n (a*) = if 

KnMI = exp [-2nr (a*)] . (47) 

Therefore, Eq. (|T9|) will hold everywhere under the barrier (ai < a < a 2 ) only if ( n (ai) is 
small enough so that 

|Cn (oi)| < exp \-2nr (a 2 )] . (48) 

The special choice £> = corresponds to the (unique) vacuum state ( n = 0. With this 
choice, the Gaussian solution will always be well-behaved, with Re S n > everywhere. 

C. Comparison: The case of if = 

If H = 0, the barrier is infinitely high, so there is no tunneling but only the recollapsing 
universe. We use this case to illustrate the underbarrier behavior of the wave function, 
because explicit calculations can be more easily done. In the H ^ case, the potential 
V (a) is well approximated by the H = expression, V (a) = a 2 — e r , when a is sufficiently 
small, a <^ H~ x . 

The Gaussian solution for H = 0, together with the condition ReS n (a) > for all a, 
gives a unique wave function, because Eq. ( |48"D forces ( — 0. The Gaussian wave function 
with C = corresponds to the vacuum state [given by Eq. (pE^)] and is illustrated in Fig. ||. 

We now consider the wave function for a squeezed state with ( = ( (ai) ^ 0. We can 
construct the exact wave function using Eqs. (p5|), (|43|) and compare it with the Gaussian 
solution. 

An explicit calculation of the exact wave function \1/ (a, Xn) is given in Appendix A. At 
a = ai it has the form of a squeezed state with a squeezing parameter £ = C(°i)- For C not 
too close to 1, the wave function is dominated by the vacuum state and nearby states with low 
occupation numbers. However, as we go under the barrier, the contribution of highly excited 
states becomes increasingly important, and at large enough a they completely dominate the 
wave function. This happens when a exceeds the following ^-dependent bound, 

a>a max (C) = ^(c^+C^). (49) 

These excited states belong to different semiclassical branches of the wave function, and 
thus the picture of a quantum field in a semiclassical background spacetime does not apply. 
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Figure 2: The vacuum wave function (( = 0) in the a — \ space with e T = 15. (Only the mode 
Xn with n = 10 is shown, and other modes \n are omitted.) Oscillatory behavior until the turning 
point ai Ri 4 is followed by an exponential decay under the barrier (a > a\). 




Figure 3: The wave function of a squeezed state (Q = 0.5) in the a — \ space, with the same 
parameters as in Fig. [2|. The underbarrier (o > a\) wave function is dominated by highly excited 
states. 
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The wave functions of the excited states exhibit short-wavelength oscillations in the x~ 
direction, which are clearly visible in Fig. |. This signals the breakdown of the Gaussian 
approximation. Now, it can be easily verified that the value a max in Eq. (|49|) coincides with 
a* defined in Eq. ([47]) as the value of a where S n (a) changes sign. This follows from fl47|) 
after substitution of the expression 



for the Euclidean conformal time for H — 0. If £ ^ then a max is finite and ReS n (a) 
becomes negative for a > a max . Only if £ = 0, the condition of Eq. (|i~9l) holds for all a > 0. 
We find that Eq. (|i~9|) is indeed the condition of consistency of the Gaussian approximation. 

We wish to stress that the constant parameter £ = £( a i) employed in the Appendix is 
not to be confused with the instantaneous squeezing parameter £ (a) defined by Eq. (]31~1) . 
The instantaneous squeezing parameter £ (a) does not necessarily reflect the true particle 
content of the quantum state. In particular, |£ (a)\ > 1 when Re S n (a) < 0, and a squeezed 
state with |£| > 1 is ill-defined and formally resembles a state with infinitely many particles, 
(N) = oo. The interpretation of Ref. 0] would suggest a "catastrophic particle production" 
due to tunneling if the barrier is wide enough so that a<i > a max . But, in fact, the model 
with a massless conformally coupled scalar field does not have any particle production. 

A physical explanation of this result has already been partially given in Ref. 0. A 
squeezed state is a superposition of all excited states, and the amplitudes of high-energy 
excited states are exponentially suppressed. On the other hand, the low-energy states con- 
tribute an exponentially suppressed amount to the wave function because of tunneling. But 
the two exponential suppression factors have a different behavior because the exponential 
suppression of low-energy states is a-dependent. Therefore, the high-energy states (which 
were already contained in the recollapsing universe) will give a dominant contribution to the 
wave function at large enough a under the barrier. 



cosh 



i 



a 



(50) 




IV. DISCUSSION 



We have examined the solutions of the Wheeler-DeWitt equation for the FRW universe 
with a conformally coupled massless scalar field. We have shown that the wave function 
under the barrier exhibits all the signs of a "catastrophic particle production" in the sense 
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of Rubakov et al., even though there is no actual particle production in this model. The 
calculations of Rubakov et al. are formally correct, but we disagree with their interpretation. 

The wave function of the universe is a superposition of different semiclassical geometries; 
each semiclassical universe contains a different quantum state of the scalar field. If the 
quantum state of the universe is a superposition of low-energy and high-energy states of 
the scalar field, then the wave function of the expanding universe will be dominated by the 
semiclassical geometries that essentially did not tunnel, rather than by the geometry of an 
inflating universe created by tunneling. 

For example, in the massless model we might consider a quantum state of the universe 
which is a superposition of the vacuum state \l/o and of the state which is the k-th 
excited state of a single mode Xn of the quantum field. A superposition of these states such 
as + a ^n) with a small amplitude a, will be close to the vacuum state to the left of the 
barrier (in the recollapsing universe). However, to the right of the barrier the contribution 
of the state will dominate the wave function, and the inflating universe will appear to 
be in an excited state with the occupation number k in the mode Xn- 

This result should be interpreted not as a production of particles during tunneling, but 
rather as an emergence of excited states that have been already present in the recollapsing 
universe and became dominant in the expanding regime after tunneling. Had these highly 
excited states not been present, the final state would have been that of an empty inflating 
universe. Thus, to investigate the creation of the universe through quantum tunneling, the 
initial state of the recollapsing universe must be chosen correctly. 

In the model with a conformally coupled massless scalar field, there is a preferred choice 
of the quantum state of the recollapsing universe which does not contain any admixture 
of excited states. This quantum state can be identified with the vacuum state. We have 
performed an explicit calculation to demonstrate that the wave function of any (non- vacuum) 
squeezed state becomes dominated by high-energy states exactly at the same region where 
the perturbative formalism of Rubakov et al. starts to break down. 

We have also shown that the Gaussian solution of the Wheeler-DeWitt equation [Eq. (^3])] 
is equivalent to a squeezed vacuum state in the formalism of Rubakov et al. The Gaussian 
solution also manifests the domination by high-energy states, in that the wave function starts 
to grow at large X- Therefore, the condition fll9| ) that the Gaussian wave function decreases 
at large x can be used to select a quantum state of the universe which is dominated by the 
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vacuum state rather than by the admixture of high-energy excited states. 

In the companion paper |[12|| . we shall extend our conclusions to the more general case 
of a massive conformally coupled field where, unlike the case of the massless field, one 
would expect some particle creation. We shall demonstrate that the quantum state of the 
recollapsing universe can be chosen to contain a sufficiently small admixture of excited states, 
and that in this case the WKB approximation is everywhere applicable and the backreaction 
of the matter excitations on the metric is negligible. We note that the under-barrier behavior 



of a massive field has been discussed by Bouhmadi-Lopez, Garay and Gonzalez-Diaz |20 



who constructed a vacuum wave function satisfying the condition of Eq. (jl9|) in the case of 
a negative cosmological constant (H 2 < 0). The paper |2(| has a significant overlap with 



our work in 12 , and we shall comment on it there in more detail. 
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Appendix A: THE UNDER-BARRIER WAVE FUNCTION OF A SQUEEZED 

STATE 

Here we compute the wave function of a squeezed state with an arbitrary squeezing 
parameter ( ^ in the case H = 0. The potential V (a) is 

V( y a) = a 2 -e r . (Al) 

For simplicity we consider the vacuum state in all modes except one particular mode \n- 
[Simultaneous squeezed states of several modes give analogous results.] Our purpose is to 
show that the wave function under the barrier is dominated by the contribution of certain 
high-energy excited states, rather than by the vacuum solution, and to find the relevant 
range of a. 

For an excited state of the mode \ n with occupation number 2k, the turning point is at 

a (2k) = yje r + 4nhk. (A2) 
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For a given a, excited states with k > k min (a) are above the barrier, where 

kmin = ~. ■ (A3) 

Ann 

The wave function for an excited state of level 2k is [cf. Eq. (|20|) ] 



*»<»•*•> = (;*) <A4) 

where the a-dependent part ip2k (a) is given by a WKB approximation, 

{ . , 1 = exp ( —4 f a . . a/V" (a) — 4nft£; da I , a > an (2fc) , 
tA^pi^ PI ftJ °f ) V U >>\ (A5) 

V^-vw C ° S (i ^ ~ V- (a) da - f ) , 0<a<a (2fc). 

The wave function ^ (a, {x/}) is a superposition of the wave functions for excited states 
of the mode Xn (we suppress the dependence on other modes xi with I ^ n), with coefficients 
given by Eq. (H) with the squeezing parameter £, 



.w-d-ifl'^-pt-g 



£j 2**! 7(45)!! ' 1 > 

For the analysis below we will need an asymptotic formula for Hermite polynomials H n (x) 
at fixed x and large n, 

H n (x) ~ v / 2(2n)t e ^-t+°(^ 1 ) cos (-^ + :ev / 2^ + Oin- 1 ' 2 )) . (A7) 

This expression can be derived by the method of steepest descent from the integral repre- 
sentation (cf. HI, Eq. 22.10.15) 

H n (a) = — / e -( 5 +-) Trft. (A8) 



Using the Stirling formula for large factorials 

n!~V2^-J (l + O^ 1 )), (A9) 



we find 



Compared with the exponential functions of k in Eq. (|A6| ) , this is a slowly changing factor. 
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Consider the contribution of the state \2k) to the wave function, as a function of k. The 
contribution of levels k > k m i n decreases with k because of the suppression factor ( k . The 
absolute value of the contribution of a level k < k m[n can be estimated as 

\C k ^2k(a,Xn)\ ~exp \-zk- \ f (a) - Anhkda] . (All) 

V n Ja (2k) J 

Here we have defined for convenience 

z = -In |C| • (A12) 

[We have used the asymptotic Eq. ( |A7|) which is justified if k is large. There is no dependence 
on Xn in the absolute value of the wave function.] For V (a) given by Eq. ( |A1|) , the integral 
in Eq. ( |A11|) can be evaluated as 

/ J a? — al da = -J a 2 — — — cosh^ 1 — . (A13) 
ia 2 v 2 a 

The function under the exponential in Eq. ( |A11|) has a maximum at 

fc = Ma - c) = i(j?x- £ -)- < A14) 

[Note that ko > only for large enough a.] Therefore the dominant contribution to the 
wave function comes from either k = (for < a < y^cosh ^-) or from k ~ k (a, C) (for 
a > y^cosh One can see this in Fig. [3|: at progressively larger values of a, the wave 
function exhibits more oscillations in the Xn direction, which corresponds to excited states 
with different values of k = k (a, Q. 

We find that the wave function \1> (a, Xn) is dominated by the contribution of excited 
states with k ~ k (a, Q when 

a > a max = ^cosh ^- = ^ (c~^ + C^ 7 ) • (A15) 



This is the expression we needed in Sec. [II C 
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